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Abstract 

The identification of pseudospin symmetry as a relativistic symmetry of the 
Dirac Hamiltonian is used to explain the structure of radial nodes occurring in 
pseudospin doublets and to illuminate the special status of nodeless intruder 
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The concept of pseudospin symmetry (T|fJ] is based on the empirical observation of quasi- 
degenerate pairs of certain normal-parity shell-model orbitals with non-relativistic quantum 
numbers 

(nr,t,3=l+\) and (n r - 1, I + 2, j = I + ^ . (1) 

Here n r , £, and j are the single-nucleon radial, orbital, and total angular momentum quan- 
tum numbers, respectively. The doublet structure, is expressed in terms of a "pseudo" orbital 
angular momentum £ = £ + 1 and "pseudo" spin, s = 1/2, which are coupled to total angular 
momentum j = £±s. For example, (n r si/2, (n r — 1)g? 3 / 2 ) will have £ = 1, (n r ps/2, — l)/s/2) 
will have £ = 2, etc. This pseudospin symmetry has been used to explain features of de- 
formed nuclei [[|, including superdeformation Q and identical bands and to establish 
an effective shell-model coupling scheme @. In view of its central role in both spherical and 
deformed nuclei, there has been an intense effort to understand the origin of this symmetry. 
The observed reduction of pseudo spin-orbit splitting in the non-relativistic single-particle 
spectra was shown to follow from nuclear relativistic mean-fields ||,[| . The pseudospin sym- 
metry itself has been shown to arise from a relativistic symmetry of a Dirac Hamiltonian in 



which the sum of the scalar, V$, and vector, Vy, potentials cancel, Vs + Vy = [pH| ) pJ| . An 
attractive scalar and repulsive vector potentials of nearly equal magnitudes, Vs ~ — Vy, is 
an inherent feature of realistic relativistic mean fields in nuclei. Calculations employing such 
fields in a variety of nuclei have confirmed the existence of an approximate pseudospin sym- 



metry in the energy spectra and wave functions [|12i --[15|| . The conditions for an approximate 



relativistic pseudospin symmetry have been further elaborated 1TB[ including the effect of 



isospin asymmetry |T1J. Implications on magnetic properties and Gamow- Teller transitions 



in nuclei |Tj| and on nucleon-nucleus scattering have been considered 



In the relativistic pseudospin scheme [ PH| , PT]| , the non-relativistic wave functions of Eq. (|TJ) 
are associated with the upper components of Dirac wave functions which are eigenstates of a 



Dirac Hamiltonian with ~ —Vy- The pseudospin s, and pseudo orbital angular momen- 
tum I, are found to be the ordinary spin, and ordinary orbital angular momentum respec- 
tively, of the lower component of the Dirac wave functions. The underlying Dirac structure 
ensures that the wave function of the upper component of the Dirac eigenfunction has a 
spherical harmonic of rank £ = £ — 1 for aligned spin: j = £ — 1/2 = £ + 1/2, and a spherical 
harmonic of rank £ + 2 = £+1 for unaligned spin: j = £+1/2 = (£ + 2) — 1/2. This explains 
the particular angular momenta defining the pseudospin doublet of Eq. (JJ). A prominent and 
uncommon feature of this doublet is that its members have different radial wave functions. 
In particular, the state with aligned spin in Eq. ([[]) has n r nodes as compared to n r — 1 
nodes of its partner with unaligned spin. One objective of the present letter is to reveal 
the so far unexplained mechanism by which the relativistic origin of pseudospin symmetry 
ensures this particular relation between radial quantum numbers of pseudospin partners. 
A related question concerns the shell- model states with aligned spin and no nodes, e.g., 
Osi/2, 0p 3 / 2 , 0d 5 / 2 , etc. For heavy nuclei such states with large j, i.e., 0g 9 / 2 , 0hu/ 2 , Ozi 3 / 2 , 
are the "intruder" abnormal-parity states which, due to the strong spin-orbit term, intrude 
from the shell above and defect to the shell below. The original formulation of pseudospin 
was oriented to a symmetry of normal-parity shell-model orbitals only. The intruder states 
were discarded from the non-relativistic pseudospin scheme. However, that is clearly not 
a satisfactory procedure if pseudospin is identified as a relativistic symmetry of the Dirac 
Hamiltonian, since now both the normal-parity states and abnormal-parity states are eigen- 
states of the same Hamiltonian. The relativistic attributes of pseudospin symmetry will 
allow us to examine the properties and special status of these intruder levels and to explain 
why these states do not have a partner which is an eigenstate. 



Properties of Dirac wave functions in a central field were considered in pO| , pT| for 
vector potentials only. The discussion below generalizes these results to the case where 



3 



both scalar, Vs(r), and vector, Vy(r), spherically symmetric potentials are present in the 
Dirac Hamiltonian. As usual, in this case the eigenstates can be written in the form 

g K (r) and f K (r) are the radial wave functions of 
the upper and lower components respectively, Yi and \ are the spherical harmonic and spin 
function which are coupled to angular momentum j projection m. The labels k = — (j + 1/2) 
and £' — £ + 1 for aligned spin j = £ + 1/2, while k = (j + 1/2) and £' = £ — 1 for unaligned 
spin j = £ — 1/2. Setting h = c = 1, the radial Dirac equations for a single nucleon of mass 
M and total energy E are 

-r- = G + A(r) F 2a 

ar r 

^L = ! F _ B {r)G (2b) 
ar r 

where G = rg K , F = rf K and 

A(r) = [E + M + V s (r) -V v (r)] (3a) 
B(r) = [E-M- V s (r) - V v (r) } . (3b) 

We consider potentials satisfying rVs, rVy — > for r — >■ 0, and Vsr, Vy — > for r — >■ oo. For 
bound state solutions — M < E < M and the radial wave functions satisfy G(0) = -F(O) = 0, 
G(oo) = F(oo) = 0, and J£°(G 2 + F 2 ) dr = 1. Specifically, we are interested in bound Dirac 
valence states for which both the binding energy (M — E) > and the total energy E > are 
positive. For relativistic mean fields relevant to nuclei, Vs is attractive and Vy is repulsive 

with typical values V s (0) 400, V v {0) ~ 350, E ~ 900, M ~ 940 MeV. Accordingly, 

in practical applications, the quantity A(r) in Eq. (3a) is positive definite, A(r) > 0, with 
A(0) ~ 1090 MeV, while the quantity B(r) in Eq. (3b) is monotonic in r and changes sign 
from a small positive value at the origin, B(0) ~ 10 MeV, to a negative value (E — M) < 
at large r. It is now straightforward to obtain from Eq. (^) the asymptotic behavior of G 
and F, 



for r — 


— > oo 


G ~ e" Ar , 


F ~ e~ Xr , 


F/G 


for r — 


-> 


G ~ r" K , 




F/G 


for r — 


-> 


G ~ r l+K , 


F ~ r K , 


G/F 



M—E 



Sr<0, («<0) (4) 
^-r>0, («>0). 
Both G and F behave as a power law for small r and exhibit an exponential falloff for large 
r. We illustrate this behavior of G and F in Fig. la,c,e for the 2p 3 / 2 {k = —2), I/5/2 = 3) 
and 0(^9/2 (k = —5) eigenstates of a Dirac Hamiltonian with scalar and vector potentials of 



Woods-Saxon form with parameters tuned to the neutron spectra of 208 Pb |L7 



To study further properties of the radial wave functions, it is convenient to introduce 
Q = r K G and T = r~ K F. Then in the open interval (0, 00), nodes of F and G coincide with 
nodes of T and Q. From Eq. (§) get 

= t 2k A(r) T (5a) 

dr 

^ = -r-2«B(r)g . (5b) 

A number of observations follow from Eq. (Bp. First, we note that it is impossible for JF 
and Q, or F and G, to vanish simultaneously at the same point because if they did, then all 
other higher-order derivatives would vanish at that point and hence the functions themselves 
would vanish everywhere. Second, we see that a node of T corresponds to an extremum of 
Q, and a node of Q corresponds to an extremum of T (since B(r) in Eq. (5b) changes sign, 
JF can have an additional extremum where B(r) = 0, which does not correspond to a node 
of Q, but this can occur only at one point since B(r) is monotonic). It follows that the nodes 
of F and G alternate; that is, between every pair of adjacent nodes of F (or G) there is one 
node of G (or F). If we let r\ be a node of T and r 2 be a node of Q, then the nature of the 
extrema at these points is determined from the second derivatives 

d 2 g 



dr 2 
d 2 T 



-A(ri)B(ri) Q(ri) where F{r\) = (6a) 



dr 2 



= -A(r 2 ) B(r 2 ).F(r2) where G(r 2 ) = . (6b) 

r=r 2 
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As bound states, both Q and T vanish at r = oo and their extrema are concave towards the 
r-axis. Therefore, the extrema at the nodes T\ or r 2 are minima (maxima) if the functions Q 
or T are negative (positive) respectively at these points. It follows from Eq. (ffl) that nodes 
of T and Q can occur only where the product A(r)B(r) > is positive. Since for practical 
applications A(r) > 0, this condition reduces to 

B(r) > at nodes of T and Q . (7) 

The combination Vg + Vy appearing in B(r), Eq. (3b), is the average potential felt by the 
nucleon in the non-relativistic limit. Therefore, the condition of Eq. ([/]) is similar to the 
statement in the non-relativistic case, that nodes of the radial wave function can occur only 
in the region of classically allowed motion, i.e., where the kinetic energy is positive. 

We now use the above results to obtain a relation between the radial nodes of F and G. 
For that purpose we consider the equation for GF = QT as derived from Eq. (fj), 

(GF)' = A(r) F 2 - B(r) G 2 . (8) 

For large r, (GF)' ~ (E + M)F 2 — (E — M)G 2 > is positive, since the binding energy 
(M - E) > for bound states. At small r, (GF)' = -B(0)G 2 < for k < 0, while 
(GF)' = A(0)F 2 > for k > 0, by employing Eq. (|). Since GF vanishes both at r = 
and r = oo we see that it is an increasing negative function at large r, while at small r, GF 
is a decreasing negative function for k < and an increasing positive function for k > 0, 

for r — ► oo GF < (9a) 
for r — ► GF < (k < 0) (9b) 
for r — > GF > (k > 0) (9c) 

consistent with Eq. (||). Furthermore, since A(r) > and by using Eq. ([7]) we find that 
(GF)' =-B(r 1 )G 2 (r l ) <0 where F(n) = 

r=n ( 10 ) 

(GF)' = A(r 2 )F 2 (r 2 ) >0 where G(r 2 ) = . 

r=r 2 



Thus GF is a decreasing function at the nodes of F, and an increasing function at the nodes 
of G. Exploiting all these derived properties, we observe that for k > 0, GF is positive 
at small r and negative at large r, and hence has an odd number of zeroes. By Eq. ( |T0D 
the first and last zeroes of GF correspond to nodes of F, and since the nodes of F and G 
alternate, then the number of nodes of F exceed by one the number of nodes of G. On the 
other hand, for k < 0, GF has the same (negative) sign near both end points, and hence has 
an even number of zeroes. By similar arguments we find that in this case the first and last 
zeroes of GF are nodes of G and F respectively, and that G and F have the same number 
of nodes. These properties of GF are illustrated in Fig. lb,d for the 2p 3 / 2 (k = —2), I/5/2 
(k = 3) states. Altogether we have, 

K < : n F = n G 

(11) 

K > : n F = n G + 1 

where np and n G denote the number of internal nodes of F and G respectively. For k > 0, 
the first and last nodes (considering F and G together) are F nodes. For k < 0, the nodes 
of G precede those of F as r increases. The same results can be obtained by considering the 
Ricatti equation for the ratio F/G and its asymptotic values, as shown in [^,0] for vector 
potentials. 

As hinted in Eq. ( p]) , the case when n F = requires a special attention. Indeed, it can 
be shown that when the radial parts of the wave functions G and F have no nodes, the 
corresponding bound states can appear only in the j = £ + 1/2 state (k < 0), but not in the 



j = £ — 1/2 state (k > 0). This result was obtained previously in [^2| for attractive scalar 
potentials of Coulomb type and separately for vector potentials of similar type. It is possible 
to generalize this result to the scalar and vector potentials under discussion. To prove the 
statement we note that according to Eq. (9a), GF is negative for r sufficiently large. If both 
G and F have no nodes, then the relation GF < holds for any r. Since by Eqs. (9b, c), 
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GF is negative near the origin only for states with k < 0, it follows from continuity that 
bound states without nodes (rip = ric = 0) must have k < 0, and hence appear only in the 
j = £ + 1/2 state, but not in the j = I — 1/2 state. An example is shown in Fig. If for the 
case of the 0g 9 / 2 (k = —5) intruder state. 

We have obtained in Eq. (|?D that the quantity B(r) is positive at the nodes of F and G. 
Even if the potentials Vg and Vy support only eigenstates with no nodes, there must still be 
a region where B(r) > 0, in order that bound states exist. This can be inferred from the fact 
that as shown GF is an increasing negative function at large r, and to enable it to vanish 
at r = 0, its derivative (GF)' must change sign from positive to negative in some region. A 
glance at Eq. (|8|) shows that since A(r) > 0, a necessary (but not sufficient) condition for 
(GF) to become negative is that 

B(r) > for some r . (12) 

The above condition, — [Vg(r) + Vy(r) ] > M — E, means that in order that bound states 
exist, there has to be a region where the depth of the average attractive single- nucleon 
potential is larger than the binding energy. 

All of the above results are relevant for understanding properties of states in the rel- 
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ativistic pseudospin scheme |l0|Jll| . The generators for the relativistic pseudospin SU(2) 
algebra, S^, which commute with the Dirac Hamiltonian, [H , S^] = 0, for the case when 
Vs — —Vy are given by 

IS.. d\ I TLS..TL n \ 

(13) 

where s M = cr^/2 are the usual spin generators, the Pauli matrices, and U p = 17 p — is 
the momentum- helicity unitary operator introduced in |§. If in addition the potentials are 
spherically symmetric, the Dirac Hamiltonian has an additional invariant SU(2) algebra, 
[H , Lfj,] = 0, with the relativistic pseudo-orbital angular momentum operators given by 
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(i 



1 ' u p i^u p o ^ 



(14) 



where = rx p are the usual orbital angular momentum operators. The sets js 1 ^, s M , s^i 



and < L 



form two triads of SU (2) algebras. The and operators act on the 



four- components Dirac wave functions. The and operators form the non-relativistic 
pseudospin and pseudo-orbital angular momentum algebras respectively, and act on the 
upper components of the Dirac wave functions. The and £^ act on the "small" lower 
components of the Dirac wave functions. The total angular momentum operators = 
+ S p have block diagonal form with entries =1^ + 5^ = U p ( £^ + ) U p = £ /jl + s^. 
The eigenfunctions of the Dirac Hamiltonian are also eigenf unctions of S • S, L • L, and 
J • J, with eigenvalues s(s + 1), £(£ + 1) and j(j + 1) respectively. The pseudospin s — 1/2 



and pseudo-orbital angular momentum £ are seen from Eqs. (p!3|)-(p^) to be the ordinary 
spin and ordinary orbital angular momentum of the lower component of the Dirac wave 
functions. This algebraic structure implies that the two states in the pseudospin doublet 
will have Dirac wave functions of the form 



^ K<0,m 



' G K<0 (r) 
r 

'Gv >0 (V 



m i 



: F K<0 (r) 
r 



Fix] 



(?) 

m 



<o, j 



(15a) 



0") ifjW 



■[Y iX ] 



+ 1>0, j" = ^+2 • (15b) 



The members of the doublet share a common pseudo orbital angular momentum £, which 
is coupled to a pseudospin s = 1/2. The state with negative k < involves unaligned 
pseudospin, j = £ — 1/2, while its partner with positive k' > involves aligned pseudospin, 

y = £ + i/2. 

In the pseudospin symmetry limit the two states in Eq. (15) form a degenerate doublet 
(S = 1/2), and are connected by the pseudospin generators of Eq. flT3p. The corre- 
sponding upper components are a doublet with respect to the set (the non-relativistic 



pseudospin algebra). Since the latter, by definition, intertwine space and spin, they 
can connect states for which the upper components have different radial wave functions, 
G K< o(r) 7^ G K i >0 (r). On the other hand, the corresponding lower components are a doublet 
with respect to the ordinary spin s^, and hence, in the pseudospin limit, their radial wave 
functions are equal up to a phase, 



K<0 



V) = F K , >0 (r) . (16) 



In particular, F K<0 (r) and F K / >0 (r) have the same number of nodes, which we denote by 
n r . If we now use the result of Eq. (0), we find for n r ^ that G K< o(r) in Eq. (15) 
has also n r radial nodes, while Cr re />o(r) has n r — 1 nodes. This explains the structure of 
nodes in the pseudospin doublets of Eq. ([!]). The simple relation in Eq. (|16|) between the 
radial wave functions of the lower components of the two states in the doublet, dictates this 
particular relation between the radial nodes of the corresponding upper components. This 
result cannot be obtained if one considers just the non-relativistic pseudospin algebra, s^, 
and is a direct outcome of the behavior of nodes of Dirac bound states and the identification 
of pseudospin as a relativistic symmetry of the Dirac Hamiltonian. 

As we have shown, bound Dirac states, for which both the upper and lower components 
have no nodes (n r = 0) can occur only for k < and not for k' > 0. From Eq. (|15a ) we 



find that such states have pseudo-orbital angular momentum I and total angular momentum 
j — £—1/2. As mentioned, these intruder states are ignored in the non-relativistic pseudospin 
scheme, and it is only the relativistic interpretation of pseudospin symmetry, combined with 
known properties of Dirac bound states, which enable a classification for these states, as 
well as provide a natural explanation why these states do not have a pseudospin partner 
which is an eigenstate of the Hamiltonian. 

The exact pseudospin limit requires that Vs(r) = — Vy{r), which implies that B(r) = 
E — M in Eq. (3b). It is clear that under such circumstances the condition of Eq. fll2|) cannot 
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be fulfilled for bound states with positive binding energy M — E > 0. This explains why 
in the exact pseudospin limit, there are no bound Dirac states and, therefore, by necessity 
the pseudospin symmetry must be broken in nuclei. Nevertheless, a variety of realistic 
mean field calculations show that the required breaking of pseudospin symmetry in nuclei is 
small |12-[T5|. Quasi-degenerate doublets of normal-parity states and abnormal-parity levels 
without a partner eigenstate persist in the spectra, and the relation of Eq. ( pf ) is obeyed to 
a good approximation, especially for doublets near the Fermi surface. As discussed, these 
features are sufficient to ensure the observed structure of nodes occurring in pseudospin 
doublets and the special status of intruder levels in nuclei. 

In summary, we have shown that identification of pseudospin as a relativistic symmetry 
of the Dirac Hamiltonian provides a natural explanation for the structure of radial nodes 
occurring in pseudospin doublets of normal-parity states. The key point in this explanation 
is that pseudospin symmetry implies a simple relation between the radial wave functions 
of the lower components of the two states in the doublet, and those in turn govern the 
radial nodes of the corresponding upper components. The intruder abnormal-parity states 
which have so far been discarded in non-relativistic treatments, can be accommodated in 
the relativistic pseudospin scheme and are assigned a pseudo-orbital quantum number £ and 
j = 1 — 1/2. General properties of Dirac bound states provide a natural explanation for 
why these nodeless states do not have a partner which is an eigenstate of the Hamiltonian. 
It is gratifying to note that characteristic features {e.g. radial and angular momentum 
quantum numbers) of states in the non-relativistic pseudospin scheme, which seem at first 
ad- hoc without an apparent reason, receive a proper justification once the relativistic origin 
of pseudospin symmetry in nuclei is taken into consideration. 
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FIGURES 



Figure 1. (a) The radial upper component (G) and lower component (F) in (Fermi) ~ 1 / 2 and 
(b) the corresponding product GF in (Fermi) -1 of the 2p 3 / 2 (ft = —2) state, (c) and (d) 
The same for the I/5/2 = 3) state, (e) and (f) The same for the 0g 9 / 2 (k, = —5) state. 
All states are eigenstates of a Dirac Hamiltonian with scalar (S) and vector (V) potentials: 
V s , v (r) = a s , v [l + exp(^)j _1 with a s = -358, a v = 292 MeV, R = 7 fm, a = 0.6 fm. 
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